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Advances  in  the  stability  and  efficiency  of  electronic  structure  codes  along  with  the  increased  performance 
of  commodity  computing  resources  has  enabled  the  automated  high-throughput  quantum  chemical  anal¬ 
ysis  of  materials  structure  libraries  containing  thousands  of  structures.  This  allows  the  computational 
screening  of  a  materials  design  space  to  identify  lead  systems  and  estimate  critical  structure-property 
limits  which  should  prove  an  invaluable  tool  in  informing  experimental  discovery  and  development 
efforts.  Here  this  approach  is  illustrated  for  lithium  ion  battery  additives.  An  additive  library  consisting  of 
7381  structures  was  generated,  based  on  fluoro-  and  alkyl-derivatized  ethylene  carbonate  (EC).  Molecu¬ 
lar  properties  (e.g.  LUMO,  EA,  fi  and  r/)  were  computed  for  each  structure  using  the  PM3  semiempirical 
method.  The  resulting  lithium  battery  additive  library  was  then  analyzed  and  screened  to  determine  the 
suitability  of  the  additives,  based  on  properties  correlated  with  performance  as  a  reductive  additive  for 
battery  electrolyte  formulations. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

For  lithium  ion  batteries,  the  electrolyte  typically  consists  of  one 
or  more  lithium  salts  dissolved  in  an  aprotic  solvent,  often  with 
at  least  one  additional  functional  additive.  Additives  are  included 
in  electrolyte  formulations  to  perform  specific  functions.  A  com¬ 
mon  function  is  to  enhance  electrode  stability  by  facilitating  the 
formation  of  the  solid/electrolyte  interface  (SEI)  layer.  Electrolyte 
additives  include  cyclic  carbonates,  such  as  vinylene  carbonate, 
vinyl  ethylene  carbonate,  fluoro  ethylene  carbonate,  and  alkyl 
carbonate  derivatives  [1-3].  Recently,  fluorinated  additives  have 
attracted  attention  in  light  of  their  ability  to  increase  the  chemical 
stability  of  electrolytes  against  high  potentials  [4].  The  discovery 
and  development  of  additives  to  increase  battery  performance  is 
of  critical  importance,  yet  such  efforts  often  involve  a  hit-or-miss 
empirical  investigative  process.  The  computational  screening  of 
chemical  design  space  for  promising  additives  has  the  potential  to 
greatly  accelerate  this  time  time-consuming  process.  Furthermore, 
theoretical  insight  into  the  effects  of  fluorination  and  or  alkylation 
on  fundamental  anode  SEI  additive  properties  would  be  useful  in 
guiding  experimental  efforts. 

The  use  of  virtual  structure  libraries  for  computational  screen¬ 
ing  to  identify  lead  systems  for  further  investigation  has  become  a 
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standard  approach  in  drug  design.  [5,6]  Transferring  this  paradigm 
to  challenges  in  material  science  is  a  recent  possibility  due  to 
advances  in  the  speed  of  computational  resources  and  the  effi¬ 
ciency  and  stability  of  materials  modeling  packages.  In  particular, 
electronic  structure  codes  are  extremely  robust  for  standard  types 
of  analyses,  usually  requiring  no  user  intervention  once  the  phys¬ 
ical  system  and  parameters  have  been  set  and  the  calculation 
initiated.  This  makes  it  possible  for  individual  calculation  steps 
to  be  executed  in  sequence  comprising  a  high-throughput  quan¬ 
tum  chemistry  workflow,  in  which  material  systems  of  varying 
structure  and  composition  are  analyzed  in  an  automated  fashion 
with  the  results  collected  in  a  growing  data  record.  This  record  can 
then  be  sorted  and  mined  to  identify  lead  candidates  and  establish 
critical  structure-property  limits  within  a  given  chemical  design 
space.  To-date,  only  a  small  number  of  studies  have  been  reported 
in  which  quantum  chemical  calculations  are  used  in  a  high- 
throughput  fashion  to  compute  properties  and  screen  for  optimal 
materials  solutions.  [7,8-1 1  ]  However,  with  time  high-throughput 
computational  screening  will  become  central  to  advanced  materi¬ 
als  research. 

In  the  current  work,  high-throughput  quantum  chemistry  and 
virtual  screening  is  demonstrated  to  explore  the  effects  of  fluorina¬ 
tion  and  alkylation  on  the  properties  of  ethylene  carbonate  (EC),  as 
the  reference.  A  virtual  structure  library  based  on  derivatized  EC  is 
enumerated  and  each  structure  is  analyzed  using  quantum  chem¬ 
ical  methods.  The  resulting  database  consisting  of  7381  structures 
is  then  analyzed  to  determine  variable  limitations  and  identify  sys- 
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Fig.  1.  Materials  discovery  and  optimization  scheme  employing  high-throughput  quantum  chemical  analyses  and  virtual  material  screening. 


terns  with  properties  characteristic  of  effective  graphite  anode  SEI 
additives. 

2.  Methods 

The  results  summarized  in  this  work  were  obtained  using  the 
VAMP  semiempirical  electronic  structure  package  [12].  Calcula¬ 
tions  were  carried  out  using  the  PM3  semiempirical  Hamiltonian 
[13].  The  sequence  of  geometry  optimization  and  single-point 
energy  calculations  for  neutral  and  charged  species  to  obtain  the 
HOMO  and  LUMO  energies,  vertical  ionization  potentials,  elec¬ 
tron  affinities,  and  chemical  hardness  was  captured  in  a  workflow 
pipeline  using  the  Materials  Component  Collection  4.4  within 
Pipeline  Pilot  (Version  7.5)  [14].  The  Pipeline  Pilot  workflow  plat¬ 
form  allowed  the  generation  and  automated  quantum  chemical 
analysis  of  an  anode  SEI  additive  structure  library  containing  7381 
stereochemically  unique  structures.  The  SEI  additive  library  was 
generated  using  an  exhaustive  combinatorial  R-group  enumeration 
scheme  (array-based  enumeration).  Both  the  electronic  structure 
code,  VAMP,  and  the  workflow  platform,  Pipeline  Pilot,  are  com¬ 
mercially  available  from  Accelrys,  Inc. 

3.  Results  and  discussion 

High-throughput  quantum  chemical  analysis  and  virtual 
screening  promises  to  revolutionize  materials  discovery  for 
advanced  applications  such  as  lithium  ion  batteries.  The  pro¬ 
cess  involves  the  generation  of  materials  structure  libraries 
(molecules,  ensembles,  surfaces),  followed  by  high-throughput 


quantum  chemical  analysis  to  compute  key  intrinsic  properties  and 
reaction  energetics  directly  related  to  battery  component  perfor¬ 
mance  and  compatibility.  The  resulting  virtual  materials  database  is 
a  unique  and  powerful  resource,  allowing  the  identification  of  opti¬ 
mal  structures/formulations  to  advise  experimental  efforts  through 
data  mining  and  screening  of  chemical  design  space  in  a  manner 
similar  to  the  scheme  outlined  in  Fig.  1.  In  the  current  work,  steps 
1  though  5  of  the  computational  materials  screening  workflow  is 
illustrated  for  lithium  ion  battery  additives. 

Semiempirical  electronic  structure  techniques  enable  the  quan¬ 
tum  mechanical  prediction  of  geometry,  electronic  structure  and 
fundamental  molecular  properties  at  a  small  fraction  of  the  cost 
of  traditional  ab  initio  methods.  The  speed,  versatility  and  estab¬ 
lished  accuracy  of  these  methods  make  them  well  suited  for 
high-throughput  quantum  chemical  analysis,  in  which  hundreds 
to  thousands  of  structures  can  be  analyzed.  In  the  present  work 
the  PM3  semiempirical  Hamiltonian  is  used  [13]  which  has  been 
shown  to  be  useful  in  predicting  the  reactivity  of  fluorinated  struc¬ 
tures  [15]  and  the  reduction  behaviour  of  electrolyte  components 
[16]. 

In  lithium  ion  batteries  utilizing  a  graphite  anode,  cyclic  car¬ 
bonates  are  often  used  as  anode  SEI  additives.  They  are  selected 
with  the  aim  of  controlling  the  chemistry  at  the  anode/electrolyte 
interface.  An  SEI  film  is  formed  during  the  first  charge  cycle.  The 
initiation  step  leading  to  anode  SEI  formation  is  electron  transfer 
to  the  SEI  forming  species  resulting  in  a  concerted  or  multi-step 
decomposition  reaction  producing  the  passivating  SEI  layer  at  the 
graphite-electrolyte  interface  [3,17,18].  Therefore,  an  important 
requirement  for  electrolyte  additives  selected  to  protect  the  elec- 


Fig.2.  Molecular  core  and  R-group  structures  used  to  enumerate  the  alkyl-  or  fluoro-derivatized  EC  structure  library:  (A)  ethylene  carbonate  (EC)  core  structure  indicating 
the  R1-R4  group  positions  and  (B)  the  various  alkyl  R-group  scaffolds  with  Z  denoting  the  connection  point  and  X  =  F  or  H. 
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Table  1 

Ethylene  carbonate  (EC)-based  additive  library  characterized  in  terms  of  relative 
property  minimum  (Min),  maximum  (Max),  mean  (Mean)  and  standard  deviation 
(Std.  dev.). 


Relative  property 

n  =  7381 

Min 

Max 

Mean 

Std.  dev. 

HOMO  (eV) 

-1.75 

0.64 

-0.56 

0.46 

LUMO  (eV) 

-3.42 

0.27 

-1.81 

0.85 

Ionization  potential  (IP)  (eV) 

-1.19 

2.20 

0.27 

0.52 

Electron  affinity  (EA)  (eV) 

-0.01 

4.13 

2.35 

0.88 

Dipole  moment  (/x)  (D) 

-4.70 

2.81 

-0.86 

1.39 

Polarizability  (or)  (A3) 

0.00 

29.00 

13.88 

4.49 

Chemical  hardness  (ij)  (eV) 

-1.69 

0.02 

-1.04 

0.32 

trolyte  solvent  from  decomposition  would  be  a  higher  reduction 
potential  than  the  solvent.  Chung  et  al.  [16]  have  shown  that 
trends  in  the  lowest  unoccupied  molecular  orbital  (LUMO)  calcu¬ 
lated  using  semiempirical  quantum  chemical  methods  reflected  the 
empirical  reduction  potentials  for  electrolyte  components.  This  is  in 
accordance  with  Koopman’s  theorem,  in  which  the  negative  value 
of  a  system’s  LUMO  energy  is  an  approximation  of  its  electron  affin¬ 
ity  [19].  Effective  anode  SEI  additives  should  possess  a  low  LUMO 
energy  or  a  correspondingly  high  electron  affinity  (EA).  An  addi¬ 
tional  favorable  characteristic  for  anode  SEI  additive  selection  is 
increased  reactivity  in  order  to  facilitate  good  SEI  formation.  One 
measure  of  the  relative  reactivity  of  a  series  of  related  structures  is 
their  chemical  hardness,  defined  by  the  following  equation  [20]: 

IPV  -  EAV 


where  IPV  and  EAV  correspond  to  the  vertical  ionization  potential 
and  electron  affinity,  respectively.  Molecules  with  a  large  IPV  and 
low  EAV  are  less  chemically  active,  i.e.,  they  are  chemically  hard; 
therefore  a  small  r\  would  be  a  favorable  feature  of  a  SEI  additive 
candidate.  Additionally,  an  effective  anode  SEI  additive  may  also 
possess  a  significant  electric  dipole  moment  (/x)  since  a  larger  value 
of  /a  leads  to  stronger  non-bonded  interactions  with  Li+  to  form 
stable  lithium  salts  as  SEI  film  components,  though  this  may  not 
always  be  a  requirement. 

To  investigate  the  substituent  effects  on  key  properties,  aris¬ 
ing  from  the  addition  of  alkyl-  or  fluoro-substituents  to  a  base 
ethylene  carbonate  (EC)  structure,  a  library  of  7381  stereochem- 
ically  unique  additive  structures  was  generated  via  an  R-group 
substitution  scheme.  The  chemical  space  is  defined  by  exhaustively 
combining  the  EC  core  structure,  shown  in  Fig.  2(A),  with  fluorine- 
and  hydrogen-terminated  alkyl  R-groups  shown  in  Fig.  2(B).  For 
each  of  the  R-groups  shown  in  Fig.  2(B),  the  number  of  carbon 
atoms  is  <3.  Quantum  chemical  analysis  of  these  structures  allows 
the  determination  of  property  limits  and  variability  of  this  focused 
chemistry  for  engineering  EC-derived  additives. 

Using  the  PM3  level  of  theory,  the  relative  molecular  proper¬ 
ties  for  each  structure  were  computed  as  the  difference  from  those 
computed  for  EC.  The  relative  properties  under  examination  are  the 
HOMO  and  LUMO  energies,  vertical  ionization  potential  (IPV),  ver¬ 
tical  electron  affinity  (EAV),  chemical  hardness  (77),  dipole  moment 
(/x)  and  mean  electronic  polarizability  (a).  The  additive  structure 
library  relative  property  ranges,  averages  and  standard  deviations 
are  presented  in  Table  1.  To  visualize  the  variation  and  distribution 
of  the  properties,  ordered  library  plots  and  histograms  for  relative 
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Fig.  3.  Ordered  plots  and  histograms  (left  and  right)  for  the  relative  HOMO  and  LUMO  energies  (top  and  bottom)  for  the  EC-based  additive  library.  Representative  additive 
structures  are  shown,  with  their  respective  property  values  indicated. 
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Fig.  4.  Ordered  plots  and  histograms  (left  and  right)  for  the  relative  ionization  potential  (IPV)  and  electron  affinity  (EAV)  (top  and  bottom)  for  the  EC-based  additive  library. 
Representative  additive  structures  are  shown  with  their  respective  property  values  indicated. 


HOMO  and  LUMO,  IPV  and  EAV,  ii  and  a ,  and  r)  are  given  in  Figs.  3-6, 
respectively.  In  each  of  the  figures,  the  left  plot  shows  the  variabil¬ 
ity  in  the  property  across  the  library  by  means  of  a  scatter  plot  (XY) 
with  the  library  structures  ordered  so  that  the  relative  property  of 
interest  is  increasing.  The  ordered  property  plots  are  annotated  to 
illustrate  representative  additive  structures  from  the  library.  The 
plot  on  the  right-side  in  each  figure  shows  the  property  distribu¬ 
tion  across  the  library  as  a  histogram  of  the  calculated  property.  An 
immediate  observation  is  that  within  the  chemical  space  defined  by 
the  EC  core  and  R-group  chemistry  shown  in  Fig.  2,  all  of  the  proper¬ 
ties  considered  here  vary  nearly  continuously  over  their  range,  with 
no  properties  showing  discontinuous  regions  of  property  values. 

Fig.  3  shows  the  variability  and  distribution  of  the  relative  HOMO 
and  LUMO  orbital  energies  with  respect  to  those  for  EC  across 
the  structure  library.  As  shown  in  the  figure  and  in  Table  1,  the 
HOMO  energy  values  range  from  -1.75  eV  below  to  0.64  eV  above 
the  HOMO  of  EC.  The  library  relative  HOMO  average  and  standard 
deviation  is  -0.56  eV  and  0.46  eV,  respectively.  The  relative  LUMO 
energies  span  from  -3.42  eV  to  0.27  eV,  with  a  mean  and  standard 
deviation  of  - 1 .8 1  eV  and  0.85  eV,  respectively.  For  reductive  anode 
SEI  additives,  a  lower  LUMO  energy  may  be  associated  with  facil¬ 
itating  SEI  film  production  due  to  easier  electron  transfer  at  the 
graphite  anode  surface.  As  seen  in  Table  1  and  Fig.  2,  the  majority 
of  the  library  structures  have  a  more  favorable  LUMO  energy  than 
EC.  Closer  inspection  of  the  nature  of  the  structures  across  the  rel¬ 
ative  LUMO  distribution  indicates  that  those  without  fluorine  have 
LUMO  energies  higher  than  that  of  EC.  Since  the  LUMO  is  only  an 
approximation  to  the  true  electron  affinity,  a  more  direct  descriptor 
for  additive  suitability  is  the  EAV. 


The  relative  vertical  ionization  potential  and  electron  affinity 
library  plots  and  histograms  are  presented  in  Fig.  4.  The  relative  IPV 
ranges  from  -1.19  eV  to  2.20  eV,  with  a  library  average  of  0.27  eV 
and  standard  deviation  of  0.52  eV.  Relative  values  of  EAV  are  pos¬ 
itive  with  a  maximum  of  4.13  eV,  mean  of  2.35  eV  and  standard 
deviation  of  0.88  eV.  This  indicates  that  all  SEI  additive  library  struc¬ 
tures  will  be  improved  electron  acceptors  over  EC.  However,  the 
effect  of  including  fluorine  in  the  molecular  structure  on  EAV  values 
is  striking.  The  minimum  and  maximum  relative  EAV  for  fluori- 
nated  structures  is  0.39  eV  and  4.13  eV,  respectively.  This  is  greatly 
improved,  compared  to  -0.01  eV  and  0.43  eV  for  purely  alkylated- 
EC  structures. 

Two  fundamental  properties  that  affect  the  intermolecular 
interactions  of  electrolyte  additives  or  co-solvents  are  the  electric 
dipole  moment  and  mean  polarizability.  The  plots  for  fi  and  a  are 
given  in  Fig.  5.  Compared  to  the  dipole  moment  of  EC,  the  additive 
structure  dipole  moments  vary  by  -4.70  to  2.81  D.  The  library  rela¬ 
tive  dipole  moment  mean  and  standard  deviation  are  -0.86  D  and 
1.39  D,  respectively.  A  sizeable  electric  dipole  results  in  a  strong 
dipole-cation  interaction  with  the  lithium  cations.  In  the  library 
2285  structures  have  dipole  moments  larger  than  that  of  EC.  The 
molecular  polarizability  is  related  to  the  condensed  phase  dielectric 
constant.  All  the  structures  have  increased  polarizability  compared 
to  EC  (up  to  +29.00 A3),  which  is  expected  since  all  the  structures 
have  a  more  extended  framework  than  the  EC  core. 

Chemical  hardness  is  related  to  the  stability  of  a  system,  reflect¬ 
ing  its  capacity  to  undergo  chemical  reactions.  From  Table  1,  the 
range  of  relative  chemical  hardness  is  1 .71  eV,  with  a  minimal  value 
of  -1.69  eV,  a  mean  value  of  -1.04  eV  and  a  standard  deviation  of 


1476 


M.D.  Halls,  K.  Tasaki  /  Journal  of  Power  Sources  195  (2010)  1472-1478 


Fig.  5.  Ordered  plots  and  histograms  (left  and  right)  for  the  relative  dipole  moment  (fi)  and  mean  polarizability  (a)  (top  and  bottom)  for  the  EC-based  additive  library. 
Representative  additive  structures  are  shown  with  their  respective  property  values  indicated. 


0.32  eV.  Fig.  6  presents  the  library  relative  chemical  hardness  plots. 
As  seen  in  the  figure,  chemical  derivatization  by  alkyl  or  fluorinated 
alkyl  groups  softens  the  system  compared  to  EC,  which  is  desirable 
for  effective  additives  as  it  should  relate  to  decreased  stability.  Here 
again  fluorinated  cyclic  carbonates  are  more  favorable  as  SEI  addi¬ 
tives,  with  a  minimum  hardness  of  - 1 .69  eV  compared  to  EC,  which 
is  0.93  eV  smaller  than  the  lowest  r\  for  alkylated  structures. 

Overall,  the  addition  of  fluorinated  substituents  to  the  ethylene 
carbonate  core  substantially  improves  the  molecular  properties 


that  correlate  with  reductive  SEI  additive  performance.  The  rela¬ 
tive  properties  are  predicted  to  vary  with  the  number  of  fluorine 
atoms  in  the  molecular  structure.  The  average  relative  LUMO 
energy,  electron  affinity  (EAV)  and  chemical  hardness  ( rj )  for  struc¬ 
tures  having  1,  3  and  5  fluorines  are  -0.1  OeV,  0.59  eV,  -0.45  eV 
(n(F)  =  1),  -0.85  eV,  1.43  eV,  -0.75  eV  (n(F)  =  3)  and  -1.51  eV, 
2.05  eV,  -1.05  eV  (n(F)  =  5),  respectively.  The  LUMO,  EAV  and  r)  for 
structures  with  n(F)  >  5  are  -2.27  eV,  2.82  eV  and  -1.18  eV,  respec¬ 
tively. 


Rel  Chemical  Hardness  (eV) 


Fig.  6.  Ordered  plot  and  histogram  (left  and  right)  for  the  relative  chemical  hardness  (r])  for  the  EC-based  additive  library.  Representative  additive  structures  are  shown  with 
their  respective  property  values  indicated. 
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Table  2 

Relative  properties  (molecular  volume,  HOMO,  LUMO,  dipole  moment  (/z),  mean  polarizability  (a),  vertical  IP  (IPV),  vertical  electron  affinity  (EAV)  and  chemical  hardness  (rj)) 
with  respect  to  ethylene  carbonate  (EC)  values  for  butyl  sultone  (BS),  propylene  carbonate  (PC),  fluoroethylene  carbonate  (FEC)  and  structure  3353  computed  at  the  PM3 
level  of  theory. 


Structure 

Formula 

MW 

Rel.  volume 

(A3) 

Rel.  HOMO 
(eV) 

Rel.  LUMO 
(eV) 

Rel.  /z  (D) 

Rel.  a  (A3 ) 

Rel.  IPV  (eV) 

Rel.  EAV  (eV) 

Rel.  rj  (eV) 

Butyl  sultone  (BS) 

c4h8so3 

138.19 

84.03 

0.68 

-1.17 

0.70 

4.10 

-0.72 

1.57 

-0.86 

Propylene  carbonate  (PC) 

c4h6o3 

102.09 

12.35 

0.09 

0.06 

0.26 

1.76 

-0.23 

0.30 

-0.26 

Fluoro-EC  (FEC) 

c3h3fo3 

106.05 

5.49 

-0.26 

-0.27 

-0.41 

0.18 

0.16 

0.42 

-0.13 

3353 

c12h13f7o3 

338.22 

138.23 

-0.17 

-2.20 

-0.46 

16.89 

-0.63 

2.74 

-1.69 

Structure  3353  (shown  in  Fig.  7)  is  a  Pareto-optimal  candidate,  simultaneously  t  (EAV  and  /z)  and  J  rj. 


For  a  given  materials  application,  finding  exemplary  structures 
involves  simultaneously  satisfying  a  number  of  target  objectives, 
in  this  case  based  on  the  fundamental  properties  related  to 
their  anode  SEI  additive  performance.  Multi-objective  solutions 
represent  a  trade-off  between  objectives,  with  one  class  being 
Pareto-optimal  solutions  [21  ].  Pareto-optimal  solutions  are  defined 
as  a  set  of  solutions  which  are  non-dominated,  such  that  it  is  not 
possible  to  improve  one  property  without  making  any  other  prop¬ 
erty  worse  [10,22]. 

As  noted  above,  for  reductive  anode  SEI  additives,  the  major 
selection  criterion  is  that  the  species  be  reduced  at  the  anode  before 
the  other  electrolyte  components  in  the  system.  Ease  of  reduction 
is  correlated  with  the  molecular  electron  affinity,  EAV.  It  follows 
that  in  selecting  for  noteworthy  SEI  structures  maximizing  EAV  is 
a  key  objective.  A  second  objective  in  selecting  SEI  additives  may 
be  to  minimize  the  chemical  hardness,  rj,  which  will  correlate  to 
a  more  reactive  species.  The  electric  dipole  of  the  system,  /x,  will 
determine  the  strength  of  the  electrostatic  interactions  between 
the  additive  and  Li+.  Screening  the  7381  EC-based  structures  for  the 
Pareto-optimal  solutions  that  maximize  EAV  and  /x,  while  simulta¬ 
neously  minimizing  rj  results  in  a  set  of  structures  that  represent 
the  best  compromise  in  terms  of  property  trade-off.  An  exemplary 
structure  from  the  Pareto-optimal  set  is  additive  3353  from  the 
library,  whose  properties  are  given  in  Table  2  and  its  molecular 
structure  is  shown  in  Fig.  7.  This  structure  has  Ri  =CEl3,  R2  =  C3El5, 
R3  =  C3F7,  R4  =  C2EI5  attached  to  the  EC  core.  The  EAV  and  rj  for  the 
structure  3358  is  computed  to  be  2.74  eV  above  and  1.69  eV  below 
those  of  the  EC  reference. 

To  gain  confidence  in  the  SEI  additive  screening  approach  pre¬ 
sented  here,  the  PM3  computed  relative  properties  for  butyl  sultone 
(BS),  propylene  carbonate  (PC)  and  fluoroethylene  carbonate  (FEC) 
are  also  presented  in  Table  2  for  comparison.  Experimentally,  butyl 
sultone  (BS)  has  been  reported  to  be  a  highly  effective  graphite 
anode  SEI  forming  additive  in  electrolyte  formulations  involving 
propylene  carbonate  (PC)  as  a  co-solvent.  [23]  The  use  of  BS  as  an 
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Fig.  7.  Molecular  structure  for  EC-based  additive  candidate  3353  representing  a 
Pareto-optimal  structure,  simultaneously  f  (EAV  and  /z)  and  4  rj. 


electrolyte  additive  overcame  performance  issues  seen  with  the 
application  of  PC  as  a  co-solvent,  such  as  the  loss  of  discharge 
capacity  and  decrease  of  cycling  stability.  BS  has  a  higher  reduction 
potential  than  PC,  thus  undergoing  reductive  decomposition  faster 
and  forming  a  passivating  SEI  film  that  suppresses  subsequent  PC 
decomposition.  A  comparison  of  the  computed  properties  for  BS 
and  PC  listed  in  Table  2  supports  the  good  performance  of  BS  as 
an  SEI  additive.  BS  is  predicted  to  be  more  easily  reduced  than  PC, 
having  an  electron  affinity  1.27  eV  larger  than  that  computed  for 
PC.  The  chemical  hardness  of  BS  is  0.60  eV  lower  than  that  of  PC, 
suggesting  it  would  be  more  reactive,  facilitating  SEI  formation. 

Fluoroethylene  carbonate  (FEC)  has  also  been  used  experimen¬ 
tally  as  an  SEI  additive  in  electrolytes  for  lithium  ion  batteries 
[24,25].  Using  an  ethylene  carbonate  (EC)  containing  solvent,  the 
addition  of  FEC  at  10%  and  30%  levels,  shifted  the  onset  of  SEI  for¬ 
mation  to  higher  potentials  by  +0.25  eV  and  +0.38  eV,  compared 
to  0.4  eV  for  the  base  solvent  vs.  Li/Li+.  The  relative  properties  for 
FEC  listed  in  Table  2  provide  a  rationalization  for  this  observation. 
As  seen  there,  FEC  has  a  higher  electron  affinity  (+0.42  eV)  and  a 
lower  chemical  hardness  (-0.13  eV)  than  EC,  suggesting  superior 
SEI  forming  behaviour. 

The  results  presented  here  introduce  a  new  and  powerful 
approach  for  exploring  the  property  limits  of  structural  motifs  for 
lithium  battery  electrolyte  additives.  High-throughput  computa¬ 
tional  screening  has  the  potential  to  dramatically  reduce  the  time 
and  effort  for  evaluating  new  synthetic  directions  for  anode  SEI 
formation  additives,  given  that  the  time  required  to  calculate  the 
key  molecular  properties  needed  for  relative  screening  was  less 
than  ca.  20  s  per  structure.  Such  screening  studies  can  help  reduce 
costs  by  informing  experimental  efforts  allowing  them  to  focus  on 
the  most  productive  candidates.  Our  findings  suggest  the  improved 
performance  for  additives  incorporating  fluorine  and  are  consistent 
with  the  limited  experimental  reports  of  fluorinated-SEI  additives 
reported  in  the  literature  [2,24,25,26]. 

The  prospect  of  automating  complex  material  modeling  work- 
flows  to  explore  formulation  and  chemical  design  space  has  the 
potential  to  rapidly  accelerate  materials  discovery,  analysis  and 
optimization.  This  approach  is  applicable  to  material  systems  rang¬ 
ing  from  discrete  molecules  and  clusters,  to  periodics  such  as 
surfaces  and  crystals.  In  addition  to  identifying  specific  lead  candi¬ 
dates  to  pursue  experimentally,  when  analyzed  using  data  mining 
and  reduction  techniques,  the  resulting  aggregate  of  data  from 
high-throughput  modeling  studies  will  reveal  structure-property 
trends  and  underlying  physical  patterns. 

4.  Conclusions 

The  use  of  high-throughput  quantum  chemistry  to  analyze  and 
screen  a  materials  structure  library  representing  a  well  defined 
chemical  design  space  has  been  applied  to  fluoro-  and  alkyl- 
derivatized  ethylene  carbonate  (EC).  The  fundamental  properties 
of  these  systems  were  evaluated  as  graphite  anode  SEI  additives 
in  lithium  ion  battery  electrolyte  formulations.  The  effect  of  flu- 
orination  leads  to  a  maximum  electron  affinity  across  the  library 
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of  4.13  eV,  compared  to  alkylation  leading  to  a  maximum  value 
of  only  0.43  eV,  relative  to  EC.  A  large  electron  affinity  is  a  crit¬ 
ical  requirement  for  improved  additive  activity.  The  library  was 
screened  to  identify  an  exemplary  Pareto-optimal  additive  struc¬ 
ture,  which  represents  a  promising  lead  for  future  investigation.  We 
propose  our  approach  as  a  new  screening  tool  in  order  to  facilitate 
the  development  of  additives  for  lithium  ion  battery  electrolytes. 
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